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f'~'~) • The production of gluons in a jet is considered in limited phase space, either with a 

cut in transverse momentum with respect to the jet axis fc^ < fc^j 11 * or with a cut in 
, absolute momentum |fc| < k cut . It is shown in the perturbative QCD calculations 

C'") ■ that in the soft limit fe<j"' — > the multiplicity distribution becomes Poissonian, 

whereas for k cut — ► the distribution remains non-Poissonian. The Poissonian 
T - H ' limit is a consequence of the soft gluon coherence in the genuine multiparticle 

^ ' correlations. We also investigate how incoherent hadronization processes could 

possibly modify the predictions for small cut-off parameters. 
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1 Introduction 

t-H . 

Multiparticle production remains an interesting testing ground for the predic- 
tions of perturbative QCD and also for the models of colour confinement. An 
' approach of great simplicity is based on the idea of duality between the prop- 

On, erties of the hadronic and partonic final state. This became first quantitative 

O 1 for the particle energy spectra which - when calculated in the modified leading 

log approximation (MLLA) - provided a good description of the datacl ( "Local 
• Parton Hadron Duality"). This poses the question whether the perturbative 

J> \ calculations are relevant to a larger class of observables including the soft par= 

1^ • ticle production. Quite a number of different observables have been calculated^ 

rS , and the comparison with experiment is encouragingly. 

d ' The standard hadronization models and the duality picture are comparable 

in the first phase of the jet evolution where secondary partons are produced 
by bremsstrahlung and pair production processes down to a relative scale of 
about i GeV. At this scale nonperturbative processes may take over as they 
are described in hadronization models. Alternatively, the perturbative phase 
may be evolved further down towards lower scales of a few 100 MeV before the 
cascade is terminated by a cut-off in transverse momentum 



k± > Qo- (1) 
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This non-perturbative cut-off represents the influence of the confinement: only 
if a gluon is radiated with larger transverse momentum it gives rise to the 
production of an extra particle. The production of hadrons may then be rather 
similar to the production of partons and in fact the data on hadron and jet 
multiplicities suggest taking the number of partons and hadrons to be equalti 

In order to test this duality picture further it is particularly important to 
investigate the final stage of the jet evolution at scales around and even below 
1 GeV. Can the decay processes involving resonance decays be described in an 
average sense by a perturbative partonic cascade? 

To approach this question it is interesting to study the peculiar aspects of 
perturbative QCD in the soft region, in particular: 

1. the running of the coupling a s (kj_) with the transverse momentum 

2. the soft gluon coherence which limits the emission angle of a soft gluon by 
the angle of the nearest colour connected parton. 

Some support for the running coupling at small scales is suggested by t 
rapid rise of multiplicity under small scale variations at small absolute scalesl 
Support for the coherence phenomenon at small scales can be drawn from 
the observation of an essentially energy independent yield of soft particles 
(below a momentum of ^200 MeV) over the full energy range explored in 
e + e~ annihilation 

In this talk I will report about a recent analysis of genuine multiparticle 
correlations in the soft limit carried out in collaboration with S. Lupia and J. 
WosiekB. The question is whether the peculiar predictions of the perturbative 
theory are followed by hadrons in this novel situation. As a consequence of 
the soft gluon coherenceQ the gluons emitted with small transverse momentum 
from the primary parton in the jet - which typically requires small angles - 
cannot radiate again. Therefore in this limit the gluon emission resembles the 
photon emission in QED in Born approximation and one expects besides the 
energy independent bremsstrahlung with a flat rapidity plateau an independent 
multiparticle emission and this corresponds to a Poissonian distribution of the 
emitted soft particles. 



2 Definition of Observables and Methods of calculation 

In order to study the multiparticle phenomena in the soft limit we introduce 
the multiplicity moments in restricted phase space 

fP>{C,Q)= [ p^(k 1 ...k q ;Q)dk 1 ...dk q , (2) 

Jt c (C,q) 



where p( q > denotes the g-particle inclusive distribution in the momenta ki at 
energy scale Q and the phase space integration is restricted by T C (C, Q) with 
a cut variable C. Two types of cuts to be applied to all particles in the final 
state have been discussed]: the momentum cut |fcj| < k cut , and the transverse 
momentum cut fcx,j < fcj"*, for i = 1, q. Clearly, the moments (|J) determine 
the multiplicity distribution of particles produced in the restricted phase space, 
and as such they provide a more differential characteristics than the global 
quantities. The latter ones are retained for the maximal k cut and k c ^ 1 at a 
given Q. 

The theoretical analysis in perturbative QCD has been carried out in two 
ways. First, analytical results have been obtained within the double logarith- 
mic approximation (DLA) M in which only the contributions of the leading 
soft and collinear singularities are kept. These results can serve to provide the 
correct qualitative picture and to test the influence of the colour coherence and 
the running coupling on the results. The DLA for the cut moments has been 
worked out completely and the results can be obtained perturbativcly to any 
order in the couplings. 

It is known since a while that the global multiplicity moments are not m= 
produced by the DLA results quantitatively. More recently it has_heen showrEJ 
that these moments are noL_well reproduced by the NLO resultsllil neither nor 
by the NNLO calculationsl!3 although a clear improvement is then obtained. 
Agreement-with the data is found finally in a numerical solution of the evolution 
equations 113. In order to obtain realistic estimates of the new cut moments we 
have therefore derived also the predictions from a Monte Carlo program (ARI- 
ADNEllJ) which is built on similar procedures than the analytical calculations. 

3 Analytical results on Multiplicity Moments 

3.1 Results in the soft limit 

The theoretical analysis starts from the equations for the generating functionals 
for the quark and gluon jets; then by the appropriate differentiation one obtains 
the equations for cumulant and factorial moments which are solved and the 
moments integrated over the regions of phase space required. 

In the case of small cut-off parameters which is of primary interest to us 
various simplifications in the analytical analysis can be applied: first, it is 
enough to take into account the terms of lowest order in the coupling, more- 
over, in this limit, the solutions for running coupling approach those for fixed 
coupling in which case close expressions are obtained. We present here only 
the main results and refer to the original paper □ for the derivation. 

Let's consider first the cylindrical phase space with small cut-off fc™*. For 
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the normalized cumulant moments wc find 



KW(X ± ,Y) = %Zl± (?±y for X X « A, Y. (3) 



in terms of the logarithmic variables 

, A = m — 

jo Qo A 



X x = hi-±-, Y = ln— , A = ln-^, (4) 



for a jet of momentum P and half opening angle 0. The numbers f q are to 
be calculated recursively (fo = 1, /i = |, /2 = §>•••)• ^ne observes that 
the moments decrease quickly with the order q corresponding to what 

has been called "linked pair ansatz" E-3 for multiparticle correlations. For the 
normalized factorial moments (F^ = 1 + =< n(n — 1) > / < n >) one 
obtains for small cut-off Xj_ 

F^(Xx,Y)^l + q -^^. (5) 

So we obtain the remarkable result that for small transverse momentum 
cut-off all factorial moments approach unity and therefore the multiplicity 
distribution becomes Poissonian. This is a consequence of the dominance of 
the single soft gluon emission at small k±, i.e., the absence of branching pro- 
cesses with secondary gluon emissions. This behaviour is just analogous to the 
usual QED bremsstrahlung and follows from_the coherence of the soft gluon 
radiation and the angular ordering conditionl3 which limits the angles of the 
secondary particle emission by the (typically small) emission angle of the first 
gluon. Namely, if this angular ordering condition is suppressed artificially the 
moments remain constant with X± and the multiplicity distribution remains 
non-Poissonian. 

On the other hand, if we carry out the same calculations for the spherically 
cut moments we obtain quite a different behaviour: the moments stay constant 
down to small cut-off parameters X = In as 

K^{X,Y) = 2*- 1 f q _ 1 /(2q-\). (6) 

Consequently, soft gluons with limited momentum k have essentially a non- 
Poissonian multiplicity distribution, while those with limited k±_ follow a Pois- 
sonian one. These calculations can be carried out for running a s as well and 
the results approach those of fixed a s in the soft limit as expected on general 
grounds. 
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3.2 General Solution in the DLA 

The multiplicity moments to arbitrary order in the coupling a s can be obtained 
by solving evolution equations in the full range of the cut-off's < X, X± < 
Y. The results reduce to the solutions for small cut-off discussed above and 
approach correctly the known results for the global moments in full phase 
space. 

The evolution equations for the unnormalized cumulant moments c^ 9 ' = 
n q in case of spherical momentum cut read 

^ h (X,Y) = I* dy^(y)[qf^- 1 \y) + f^(y)](X-y) 
Jo 

[■X rY-X+x 

+ dx dy^(y)f^ h (x,y) (7) 

JO Jx 

where f^(Y) refer to the known global unnormalized factorial moments and 
7q = 6a s /ir. From this equation the cut moments can be calculated recursively 
to any precision. 

The moments in cylindrical phase space can be obtained from those in 
spherical phase space using the equation 

c^ l {X^Y)=c^{X 1 _)+ 

(Y x ± ) [ X± 7o 2 (y)[9/ ( ^ 1} GO + f {q) (y)]dy + 

Jo 

f (Y-yh 2 (y)f^ h (X^y)dy. (8) 
Jx ± 

The results from this calculation up to 0(af ) are shown in Fig. [j] in the 
full range of cut-off parameters X, X± . The figure displays the very different 
behaviour of the moments in spherical and cylindrical phase space with the 
latter ones approaching unity for X±_ — > corresponding to the Poisson limit. 

4 Monte Carlo results 

We complement our analytical DLA calculations by a Monte Carlo computa- 
tion of the multiplicity moments. Our aim is to clarify whether the qualitative 
features of our DLA calculations are verified and to obtain quantitative pre- 
dictions at the parton level. We have chosen the ARIADNE programs which 
has a similar procedure in evolving the parton cascade and in the final cutoff 
(uh which can be chosen independently of the QCD scale A. By comparing the 
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Figure 1: DLA predictions for the cut-off dependence of the first three normalized moments 
F^ , q = 2, 3,4 for Y = 5.7. Both families, i.e. spherically (sph) and cyllindrically 
(cyl) cut moments, are shown. They coincide in the global limit (X or X± = Y), but 
have distinctly different threshold behaviour. Different lines which describe one moment 
correspond to different orders k max of the perturbation theory included k max = q, q+1, 8. 



Monte Carlo results with the data on the total multiplicities and global mo- 
ments in e + e~ annihilation we adjust the parameters of the program Qo = 0.2 
GeV and A = 0.015. With these parameters we obtain the parton level predic- 
tions for the cylindrical moments in Fig. ^|a. The moments clearly show the 
decrease towards unity for decreasing cut parameter fc™* in qualitative agree- 
ment with our analytical calculation. The absolute size, however, is lower in 
the Monte Carlo than in the analytical calculation which neglect various fi- 
nite energy effects. The moments for the spherical cut (not shown) keep some 
larger non-zero values at small cut-off. We take these findings as confirmation 
of our main analytical results on the very different behaviour of the two fami- 
lies of moments and of the approach to the Poissonian limit of the cylindrical 
moments. 

Furthermore, we studied the effect of hadronizationxin the considered ob- 
servables. The idea of Local Parton Hadron Duality El has been originally 
proposed for the single inclusive distributions. The study of the moments 
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Figure 2: (a) cut moments of order 2 (diamonds), 3 (crosses) and 4 (squares) in one hemi- 
sphere defined through the thrust axis as a function of fc^ 11 * as predicted by ARIADNE at 
parton level with parameters readjusted; (b) same as in (a), but at hadron level with default 
values of the parameters and string fragmentation. 
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considered here allows for a novel test of these ideas in the case of genuine 
multiparticle correlations. In Fig. ||b we show the results of the ARIADNE 
program with string hadronization using the standard parameters. There is the 
remarkable prediction that the moments first decrease with decreasing fc™* as 
in the parton level calculation but for fc™* < 0.8 GeV the moments rise again 
and do not approach the Poissonian limit. 

5 Conclusions and Outlook to Other Processes 

We have shown that perturbative QCD predicts the different behaviour of 
multiplicity moments in cylindrical and spherical phase space regions. The 
approach to a Poissonian in the soft limit in case of the cylindrical cut is a 
consequence of the soft gluon coherence. This behaviour is also supported by 
the Monte Carlo calculation which takes into account the most important non- 
leading effects. The Poissonian multiplicity corresponds to the independent 
emission of soft gluons distributed with a flat rapidity plateau in close similarity 
to the QED Bremsstrahlung. On the other hand, hadronization effects may 
distort the predicted behaviour. It will therefore be very interesting to study 
these moments experimentally and to find out whether, or, to what extent 
the perturbative predictions survive the hadronization process in this case of 
genuine multiparticle correlations. 

Our calculations apply to quark and gluon jets in hard collisions. If the 
(untriggered) high energy pp collisions proceed with semihard gluon exchange 
and subsequent gluon bremsstrahlung a similar behaviour is expected and the 
moments with cylindrical cut should approach a Poissonian distribution. On 
the other hand, in collisions of heavy nuclei, if there is a thermalization process 
with quark gluon plasma formation, the angular ordering in the cascade process 
is destroyed and no Poissonian limit is obtained. In this way the multiplicity 
moments could provide a new probe of thermalization. 
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